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INTRODUCTION 


A  conventional  wick  employs  surface  tension  to  raise 
fluids  against  gravity  and  is  used  in  applications  where 
only  smell  flow  rates  are  required.   The  wick  is  extremely 
simole  and  reliable  and  requires  no  external  source  of 
power.   An  example  of  current  interest  where  a  conven- 
tional wick  is  applied  is  in  the  liquid  return  path  of  a 
heat  pipe.   Heat  pipes  are  used  in  many  applications  where 
large  quantities  of  heat  must  be  removed  from  a  confined 
space.   Examples  include  heat  transfer  from  the  core  of 
a  nuclear  reactor,  the  cooling  of  large  electronic  tubes 
and  semiconductor  components,  and  heat  removal  from  trans- 
formers.  Heat  is  picked  up  by  the  working  fluid  of  the  heat 
pipe  as  latent  heat  of  vaporization  when  the  liquid  evapo- 
rates and  is  convected  along  with  the  vapor  until  the  liquid 
condenses  in  the  condensor  section  end  the  heat  is  rejected 
to  some  heat  sink.   The  liquid  must  then  be  returned  to  the 
evaporator  in  order  that  the  cycle  be  continuously  repeated 
and  this  is  usually  accomplished  with  a  conventional  wick. 
Using  this  orocess,  heat  fluxes  ctin   be  achieved  which  are 
much  higher  than  those  possible  by  conduction  through  a 
metallic  conductor  with  the  same  temperature  difference 
between  the  heat  source  and  sink..  The  limitations  uoon 
heat  transfer  rate  which  can  be  achieved  with  a  heat  pipe 


fell  Into  two  areas:  first,  the  flow  rate  of  the  liquid  re- 
turning from  the  condensor  to  the  evaporator  is  limited  if  a 
conventional  wick  is  used  as  the  return  mechanism;   second, 
the  heat  flux  in  the  evaporator  is  limited?  because  of  the 
vapor  film  which  forms  at  the  liquid-metal  interface  and  acts 
as  an  Insulator  to  increased  heat  flux* 

Limitations  on  flow  rate  In  the  conventional  wick 
occur  because  of  the  competing  effects  of  surface  tension 
and  viscous  drag.   Liquid  will  rise  through  the  pores  of  the 
wick  until  the  surface  tension  force  at  the  liquid  surface 
balances  the  gravltlonal  body  force  on  the  liquid  in  the  pore. 
To  increase  the  potential  heighth  to  which  the  liquid  can  be 
raised,  the  average  pore  diameter  of  the  wick  must  be  de- 
creased to  increase  the  ratio  of  liquid-solid  interface  to 
cross  sectional  area.   As  the  fluid  begins  to  flow  the  sur- 
face tension  forces  must  compete  with  another  pressure  loss 
due  to  viscous  drag.   Viscous  drag  increases  as  the  diameter 
of  the  pores  are  decreased.   Hence  for  a  given  heighth  to 
which  the  fluid  must  be  raised,  there  Is  a  limited  flow  rate 
which  can  be  realized. 

An  increase  in  liquid  flow  rate  by  using  an  independent 
pumping  mechanism  could  provide  an  increase  in  the  heat 
transfer  capability  of  a  heat  pipe.   However  one   of  the 
attractions  of  the  conventional  heat  pipe  is  that  the  wick 
Dumping  does  not  require  an  independent  energy  source  and  is 
highly  reliable.   But  if  the  energy  source  used  for  pumping 
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could  also  be  used  to  help  overcome  the  other  limitstiojjl  on 
the  heat  pipe  performance,  i.e.  the  maximum  heat  flux  due  to 
film  boiling  in  the  evaporator  and  condensation  in  the  cion- 
densor,  improvement  of  the  total  performance  of  the  heat  pipe 
becomes  an  attractive  possibility. 

There  has  been  investigation  of  improvement  in  evapora- 
tion and  condensation  heat  flux  by  using  electric  fields  to 
cause  the  vapor  film  to  become  unstable  and  breaK  av,ey  from 
the  metal  boiler  or  condensor  surface.  This  requires  a  high 
voltage  source  and  suggests  that  if  such  a  scheme  could  be 
applied  to  the  heat  pipe,  a  natural  area  of  investigation 
toward  improvement  of  the  total  heat  pipe  performance  would 
be  the  use  of  electric  fields  to  provide  the  pumping  force 
for  the  v:  or  king  liquid. 

The  use  of  electric  fields  in  pumping  of  liquids  can  take 
advantage  of  either  dielectrophoretlc  force  where  the  force 
is  perpendicular  to  the  electric  field  direction,  or   the  body 
force  on   the  fluid  due  to  the  ion  drag  phenomenon  where  the 
force  is  in  the  direction  of  the  applied  field.   In  either 
case  the  fluid  must  be  a  grod  insulator  and  should  resist 
electrical  breakdovn  as  well  as  possible  because  the  forces 
exerted  increase  as  the  field  increases.  A   fluid  satisfying 
these  requirements  is  freon  113  which  is  used  widely  as  a 
heat  transfer  agent. 

The  requirements  of  the  pumping  mechanism  for  the  heat 


pipe  application  are  these: 

(1)  It  must  be  possible  to  achieve  high  flow  rates  re- 
lative to  those  possible  with  a  conventional  v  ick0 

(2)  It  must  be  possible  to  raise  the  liquid  to  an 
arbitrary  helghth  in  a  gravitational  field  or  to 
move  the  liquid  against  adverse  accelerations  such 
as  those  encountered  in  a  space  application. 

(3)  The  pumping  mechanism  must  be  highly  reliable. 
(b)      Simplicity  and  ease  of  fabrication  are  necessary 

if  it  is  to  compete  with  a  conventional  wick. 

The  use  of  dielectrophoretic  force  for  pumping  fails  on 
the  second  requirement  above.   The  maximum  pressure  which  can 
be  generated  at  the  interface  of  a  liquid  is  j((-^)f2  where 
(  is  the  permittivity  of  the  liquid  and  C0  is  the  permit- 
tivity of  the  soace  adjacent  to  the  liquid  surface.   For 
freon  this  pressure  represents  less  then  6  centimeters  head 
for  the  maximum  field  which  can  be  tolerated  without  elect- 
rical breakdown.   In  addition  no  convenient  method  of  stag- 
ing the  pumps  is  available;  hence  the  pressure  rise  achiev- 
ed in  separate  stages  cannot  be  added.   Thus  the  maximum 
vertical  dimension  of  a  heat  pipe  using  this  method  of  pump- 
ing would  be  less  than  6  centimeters  if  the  condensor  were 
located  below  the  evaporator.   Such  a  limitation  would  not 
be  imposed  in  a  space  application  but  the  adverse  accelera- 
tions under  which  the  pipe  could  operate  would  be  limited. 

The  remaining  Dosslbility  for  application  of  electrical 


pumping  is  that  of  the  ion  dreg  phenomenon.   Investigation 
of  this  phenomenon  has  been  made  most  recently  by  Stuetzer 
(1,2,3)  and  Pickard.(5)   a  practical  pump  has  been  suggest- 
ed and  patented  by  the  former  and  has  been  used  to  achieve 
high  pressure  rised  in  the  fluids  Investigated.   However  the 
geometries  used  are  quite  difficult  to  fabricate  especially 
in  the  numbers  that  would  be  required  for  the  heat  pipe  wick 
apolication.   In  fact  the  requirement  of  this  application  is 
to  achieve  high  flow  rates;  pressures  need  be  only  high 
enough  to  lift  the  liquid  to  the  next  stage  of  the  pump. 
Therefore  a  simple  geometry  for  ion  drag  pumping  which  per- 
mits high  flow  rates  and  adquate  pressures  is  the  basic  re- 
quirement of  this  application.   Proposal  and  investigation 
of  such  a  geometry  is  the  subject  of  this  thesis. 

The  ion  drag  phenomenon  can  be  characterized  as  follows  : 
A  body  force  can  be  exerted  on  an  insulating  fluid  if  charged 
particles  in  the  fluid  are  pulled  along  by  an  electric  field. 
The  particles  transfer  momentum  to  the  fluid  through  friction 
as  they  move  under  the  influence  of  the  electric  field,  and 
a  pressure  gradient  is  generated,  in  the  fluid  which  results 
in  a  pumping  notion  if  the  fluid  is  permitted  to  move.   The 
magnitude  of  the  pressure  which  can  be  realized  is  limited 
by  the  electrical  breakdown  strength  of  the  fluid.   A  typi- 
cal arrangement  is  one  in  which  two  electrodes  are  used,  one 
being  a  point  at  which  corona  discharge  produces  the  nec- 
essary charged  particles  and  the  other  electrode  being  more 
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blunt  so  that  few  charged  particles  are  produced  there. 
When  a  sufficiently  high  voltage  is  placed  between  the  elect- 
rodes, a  pressure  gradient  is  generated  and  pumping  of  the 
fluid  will  usually  occur  from  the  point  to  the  blunt  elect- 
rode if  the  flow  geometry  permits.   In  order  to  produce 
pressures  large  enough  to  be  of  interest,  it'  is  necessary  to 
provide  lateral  boundries  of  insulating  material  between  the 
electrodes  to  insure  that  most  of  the  charged  particles  flow 
in  the  desired  direction  of  pumping. 

Because  of  the  requirement  for  lateral  boundries  in  the 
ion  drag  pump  the  physical  size  of  a  single  pump  is  usually 
quite  small  and  therefore  the  flow  rate  which  can  be  real- 
ized is  not  large.   In  addition,  the  physical  geometry  of 
the  pump  provides  a  large  hydrodynamic  resistance  factor 
which  produces  a  pressure  loss  depending  upon  the  flow  rate. 
As   the  flow  rate  increases  the  point  at  which  the  pressure 
loss  exceeds  the  generated  pressure  is  soon  reached. 

In  order  to  make  practical  use  of  ion  drag  pumps  for  the 
wick  application  to  a  heat  pipe,  it  is  necessary  that  one  be 
able  to  fabricate  large  numbers  of  these  pumps  in  parallel  to 
achieve  the  necessary  flow  rates.   Also  since  we  must  sacri- 
fice a  large  percentage  of  generated  pressure  to  hydrodynamic 
losses  in  the  pump  itself,  an  arrangement  for  placing  the 
paralleled  arrays  of  pumps  in  series  is  necessary.   Since  the 
objective  is  to  use  the  entire  arrangement  as  a  wick,  the  re- 
sultant pressure  from  each  stage  need  be  only  slightly  great- 


er  than  necessary  to  raise  the  fluid  vertically  through  the 
pump  dimension  and  to  the  next  stage.   Therefore  a  geometry 
will  be  sought  vihich  maxumizes  flow  rate  while  maintaining 
sufficient  pressure  to  operate  as  a  wick. 

Finally  the  objectives  of  provision  for  parallel  and 
series  operation  must  be  achieved  with  a  geometry  which  can 
be  easily  fabricated  and  whose  dimensions  may  be  easily  ad- 
justed to  achieve  the  desired  tradeoff  between  pressure  and 
flow  rate. 

The  possible  application  of  an  electrohydrodynaraic  wick 
satisfying  the  given  requirements  is  not  limited  to  heat 
pipes.   However  it  was  this  application  which  provided  the 
motivation  for  study.   Another  possible  application  would  be 
to  provide  motion  of  the  oil  inside  transformers  thus  elim- 
inating hot  spots.   Tom  Jones  has  used  the  configuration  pre- 
sented in  this  thesis  in  a  wall-less  pipe.'' ' 

The  ion  drag  pumping  of  liquids  has  been  observed  for 
many  years  but  little  application  has  been  made  of  the  phe- 
nomonen.   V.rith  the  wick  application  to  heat  pipes  in  mind, 
the  goal  of  this  investigation  has  been  to  develop  a  con- 
figuration so  simple  it  could  be  made  and  used  by  anyone 
having  requirement  for  a  similar  pumping  device,  and  to 
determine  the  characteristics  of  the  device  sufficiently 
well  that  a  systematic  design  can  be  made  for  any  required 
application. 


CHAPTER    I 
PUMP  GEOMETRIES 


Pump  Geometries  Invest igated 

The  geometry  investigated  in  most  experiments  reported 

( 1  ) 

by  Stuetzer v w  was  that  of  a  nolnt  to  ring  as  shown  in  figure 

1-1.   Corona  discharge  at  the  point  electrode  produces  free 
charge  which  is  pulled  toward  the  ring  electrode  and  the 
fluid  is  pumped  in  the  direction  shown.   High  pressure  rises 
in  the  fluid  can  be  realized  across  the  pump.   The  disadvan- 
tage of  this  geometry  is  that  fabrication  of  many  of  these 
pumps  is  difficult  because  point  electrodes  must  be  placed 
within  the  channel  at  a  precise  distance  from  the  ring.   In 
addition  the  point  electrode  and  its  support  provides  an 
additional  hydrodynamic  resistance  to  the  flow.   This  is 
significant  for  the  present  application  because  the  object 
is  to  operate  at  low  pressure  rise  and  high  flow  rates. 

The  geometry  used  by  Pickard  in  his  investigation  was 
that  of  two  grids  of  thin  wire  across  the  channel  as  shown 
in  figure  1-2.   The  same  effect  can  be  achieved  by  using  wire 
screens.   These  offer  hydrodynamic  resistance  to  the  flow;  in 

addition,  wire  or  screen  fine  enough  to  act  as  a  corona  emit- 

i 
ter  is  easily  broken  and  makes  fabrication  difficult. 

Cther  geometries  were  investigated  including  the 

straight  wire  to  ring  configuration  shown  in  figure  I-3» 

Here  the  charges  are  produced  near  the  thin  wire  and  the 

fluid  is  pumped  toward  the  ring  electrode.   This  geometry  is 

easier  to  fabricate  than  the  point  to  ring  since  the  same 

wire  can  be  placed  across  the  channels  of  adjacent  pumps  in 
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parallel.   However  the  electric  field  configuration  is  not 
symetric  about  the  channel  axis  being  stronger  at  the  sides 
near  the  wire  and  weaker  90  degrees  from  this  point  around 
the  channel.   Circulation  of  the  fluid  between  the  electrodes 
occurs  before  the  pressure  developed  is  as  high  as  desired. 
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Figure  1-2 
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Figure  1-3 
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Proposed  Pump  Geometry 

The  pump  geometry  which  was  selected  and  tested  is  as 
shown  in  figure  1-4.   It  consists  simply  of  two  thin  metal 
plates  with  a  sheet  of  plexiglr.s  sandwiched  between.   To  make 
the  pump,  a  hole  is  drilled  through  the  metal  plexiglas 
sandwich.   As  the  drill  bit  emerges  from  the  first  brass 
sheet  and  enters  the  plexiglas,  a  burr  or  projecting  sharp 
edge  is  formed.   In  freon,  testing  has  shown  that  this  is 
sufficient  to  produce  the  necessary  free  charge  as  ions 
within  the  liquid.   The  channel  through  the  plexiglas  dir- 
ects the  ions  such  that  their  motion  toward  the  other  elect- 
rode causes  a  pressure  to  be  built  up  between  points  (A)  and 
(B)  in  figure  1-4.   The  electrodes  are  secured  to  the  plexi- 
glas before  the  hole  is  drilled  so  that  the  pump  need  not  be 
realligned  after  drilling.   The  extreme  simplicity  of  the 
geometry  makes  it  easy  to  produce  pumps  with  a  variety  of 
electrode  spacings  (by  changing  the  hole  size).   By  doing  so 
it  is  possible  to  determine  optimum  dimensions  to  produce  the 
best  performance  with  respect  to  pressure  and  flow  rate. 

Parallel  operation  of  pumps  is  easily  produced  by  drill- 
ing several  holes  as  shown  in  figure  1-5 «   The  possible 
interaction  of  pumps  which  might  degrade  performance  was  in- 
vestigated to  determine  if  a  minimum  spacing  of  the  pumps  is 
required.   It  was  also  necessary  to  determine  if  slight 
differences  In  electrode  edges  would  degrade  the  overall  per- 
formance of  the  entire  array. 
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The  series  arrangement  Is  shown  in  figure  1-6.   It  con- 
sists of  several  layers  of  the  parallel  arrays.   The  regions 
between  the  arrays  are  reservoirs  for  the  liquid.  On   either 
side  of  these  reservoirs  are  exposed  electrodes  of  opposite 
polarity  and  placed  such  that  they  would  cause  a  loss  in 
pressure  if  ions  v;ere  pulled  across  this  space  in  opposition 
to  the  desired  flow  direction.  However,  since  there  are  con- 
stricting boundries  as  there  are  in  the  -Dump  region,  little 
loss  was  observed.   In  any  case,  such  loss  can  be  decreased 
by  Increasing  the  spacing  between  pump  arrays. 

One  difficult ity  with  the  given  pump  geometry  is  that 
it  is  not  easy  to  analyze  mathematically.   It  does  not  fall 
exactly  into  any  of  the  three  categories  analyzed  by  Stuetzer, 

I.e.  plane,  cycllndrical,  or  spherical.   The  emitter 
electrode  is  not  a  plane  surface,  a  point,  or  a  straight 
line  but  a  ring  at  one  end  of  the  cyclinder.   Therefore  it 
is  necessary  to  base  the  analytical  model  on  an  approximate 
geometry.   Experimental  study  of  the  actual  pump  was  used  to 
determine  what  adjustments  were  necessary  to  the  mathematical 

model  in  order  to  have  reasonable  predictions  of  performance 

i 

to  be  used    in  design  of  practical  wick  arrays    for  a   s-oeclfic 
amplication. 
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CHAPTER  II 
STATIC  PRESSURE  GENERATION 
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Stet 1c  Pressure  Generation 

The  static  model  for  the  ion  drag  device  relates  the 
external  fluid  pressure  generated  with  the  applied  voltage 
since  the  device  will  be  driven  by  a  constant  voltage  d.c. 
source.   The  geometry  explained  in  chapter  -I  is  complicated 
to  model  analytically.   Therefore  development  of  the  model 
will  be  done  as  follows.   Experiments  have  shown  that  for  a 
particular  set  of  pump  dimensions  the  pressure-voltage 
dependence  can  be  modeled  by  a  space  charge  limited  current 
and  the  problem  can  be  solved  if  a  simple  geometry  is  chosen 
as  done  by  Stuetzer' -^ ).   This  theory  is  then  modified  to 
account  for  the  observed  performance  of  the  proposed  pump 
geometry.   Experiment  has  also  shown  that  there  is  a  depend- 
ence of  generated  static  pressure  upon  the  transverse 
dimensions  of  the  nump..  This  dependence  was  investigated 
experimentally  to  determine  optimum  dimensions  for  wick 
anpllcation. 
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One  Dimensional  Static  Hod el  j 


60  \£'  Flu,  4      \£'(b) 


-o  *      x 


A   one  dimensional  model  as  shown 
in  figure  II-l  will  be  assumed.   The 
charged  particles  enter  the  fluid  at        Figure  II-l 
the  electrode  E  and  are  removed  from 

the  fluid  at  electrode  E?.   Therefore  the  charge  relaxation 
time  of  the  liquid  must  be  greater  than  the  transient  time 
of  the  particles  between  electrodes.   The  electric  field 
seen  by  the  uarticles  will  be  the  macroscopic  field  of  the 
system  end    individual  particles  will  not  influence  one 
another.   The  liquid  will  appear  as  a  continuum  to  the 
moving  particles..  It  is  assumed  that  only  one  type  of  charged 
particle  Is  present.   This  can  be  realized  if  the  electrodes 
are  not  symmetrical  or  if  the  fluid  exhibits  a  strong  tend- 
ency to  produce  only  one  type  ion  in  a  corona  discharge. 
The  electrical  force  on  a  particle  will  be 

where  q  is  the  particle  charge  and  E  is  the  local  field  at 

the  particle.   The  viscous  force  exerted  on  the  particle  by 

i 
the  fluid  is  pr ooortatlonal  to  the  velocity  of  the  particle 

relative  to  the  stationary  fluid 

The  particle  velocity  can  be  expressed  in  terms  of  the 
electric  field  since  the  viscous  enc   electrior.l  forces  must 
balance  in  steady  motion 
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The  quantity  b  Is  called  the  mobility  and  v.  ill  be  assumed 
independent  of  the  field.  If  the  fluid  is  stationary  and 
charge  diffusion  is  neglected  the  current  then  can  be  ex- 
pressed  as 

$  -    f  v   =     f  b£  '  (1) 

v;here   f  is   the   charge  density  and  a   is   the   current  density. 
The   total  body    force   exerted   on  the    fluid  per   unit   volume    is 
the    total   electric    force 

P  =    ?t 
Since    the    fluid    is    not    in  motion  there    are   no    inert ial   or 
viscous    forces   and   the   body   force    must   be   balanced  by   a 
pressure   gradient. 

P*  ~     V  p  (2) 

Surface   charges   on  the    circumference   of   the   channel  villi  be 

neglected.      Hence    the    charge    density   and   the   electric    field 

will   be   uniform  across   the    channel   and   the   electric    field 

will  have   only  an  x   component.      Gauss's    Law    is   then  written 

sUL   _  JL  .       (3) 

and  the  electric  field, is  related  to  the  voltage  through 

■cHc    ~       L  (M 

Substitute    (3)    into    (1) 

and    integrate    from   X=  o   to  x   snd  £=  /e    to    £ 
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£»-  (^    +  ^   ^ 


From  eq0   (3)  p  d£ 


(5) 


9/h (6) 

To  find  p  integrate  from  #=o,  .xand  p -' o     /••-  p 


To  introduce  voltage  Integrate  eq.  (5)  using  eq.  (M 

Z*  =  r  =   Ta  f  £e  7i 


U  =      LA  ( i_2*  ,  ,.  i^ 


*  -     <*    (-f^  ,  £e;* 


Specific  Emitter  Assumption 

If  it  is  now  assumed  that  Ee   =  o   .   This  assumption  will 
be  discussed  later.  Then 

e  -  ,f  «  (£/  (8) 

Equation  (8)  shows  that  the  pressure  difference  be- 
tween points  (A)  and  (B)  in  figure  TI-1  is  independent  of 
the  lateral  dimension  of  the  channel.   This  is  because  the 
space  charge  is  uniformly  distributed  across  the  channel  and 
there  are  no  field  components  perpendicular  to  the  x  direct- 
ion . 


?0 


Comparison  of  equation  (8)  to  the  result  obtained  by 
Stuetzer  (1)  for  a  spherical  (point  to  ring)  geometry  shows 
that  for  a  fixed  set  of  pump  dimensions,  the  relationships 
differ  by  only  a  constant.   Since  an  exact  relationship 
cannot  be  obtained  for  the  geometry  proposed  in  chapter  I, 
a  relationship  of  the  following  form  is  assumed; 

where  the  value  of  C  is  to  be  determined  by  experiment. 
In  terms  of  pressure  head       /  __  C  6 _  (  U '  V2, 

«  W        (9, 

where  h  is  the  head  developed  across  the  pump.   Here  the 

voltage  U'  is  the  difference  between  the  applied  voltage 
and  a  "cut  in  voltage"  usually  taken  as  the  minimum 
voltage  necessary  to  produce  unipolar  ionic  conduction. 

The  dependence  of  pressure  proportat ional  to  voltage 
squared  is  a  result  of  the  boundry  condition  assumption  that 
the  electric  field  at  the  emitter  is  zero.   Figure  II- 2 
shows  a  plot  of  pressure  head  developed  across  a  pump  as 
a  function  of  U  .   This  linear  relationship  is  typical 
for  most  Dumps  investigated  in  the  region  of  voltage  where 
the  pumps  are  to  be  operated.   Breaks  have  been  observed 
in  these  plots  but  these  will  be  discussed  more  in  chapter 
V.  I 

The  value  of  C  in  equation  (9)  is  determined  from  a 
plot  such  as  figure  II-2  from  the  equation 


?1 


C  ~  7.*y     x*   %hjx  do) 


44 

where  x  is  the  electrode  speclng  in  millimeters,  "3JP-  is  the 
slope  of  the  pressure-voltsge  squared  plot  in 
centimeters/( kilovolt )  . 

The  value  of  C  depends  upon  the  diameter  of  the  hole 
and  the  size  of  burr  produced  in  drilling.   However,  it  has 
been  found  that  for  most  pumps  with  spacers  of  1/16  inch 
and  holes  of  diameter  approximately  1.^0  millimeters,  the 
value  of  C  is  in  the  range  0.8  to  1.1  . 

The  spacing  of  1/16  inch  was  found  to  work  quite  well 
at  reasonable  voltages  of  10-15  KV.   Smaller  spacings  cause 
electrical  breakdown  at  voltages  only  slightly  above  the 
"cut  in  voltage"  and  thicker  spacers  required  extremely  large 
voltages  (25-30  KV)  to  give  significant  pressures.   At  the 
spacing  of  1/16  inch  breakdown  would  usually  occur  at  16  to 
19  KV  depending  upon  the  purity  of  the  freon. 
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Pressure  Heed  vs 
Voltage  squared 


Voltage  squared  (KV) 


Figure  TI-2 
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;xperi mental  Effects  of  Lateral  Pump  Dimensions 


The  preceeding  development  does  not  account  for-  any 
pressure  dependence  upon  the  lateral  dimension  of  a  pump, 
i.e.  the  diameter  of  the  hole  for  the  geometry  as  proposed 
in  chapter  I.   As  has  been  stated  the  exact  theoretical  solu- 
tion for  this  geometry  has  not  been  worked  out  so  an  experi- 
mental investigation  was  carried  out  for  a  large  range  of 
hole  sizes  with  pump  length  fixed  at  1/16  inch.   Figure  TI-3 
shows  a  plot  of  pressure  vs.  pumo  diameter  at  three  different 
voltages  for  constant  pump  length.   The  freon  used,  in  these 
experiments  was  doped  with  copper  napthenate. 

For  very  large  holes  the  generated  pressure  decreased 
with  increasing  diameter.   Tons  are  being  injected  only  at 
the  outside  edges  of  the  channel  and  it  would  be  expected 
that  for  very  large  diameters  the  pressure  due  to  drag  of 
the  ions  toward  the  collector  would  be  small. 

The  interesting  result  is  shown  by  the  jump  in  pressure 
at  a  diameter  of  about  1.9  millimeters.   For  holes  smaller 
than  approximately  1  mm.  the  pressure  does  not  change  much 
in  these  experiments.   It  is  not  practical  to  make  holes 
smaller  than  about  0.2  mm  with  drill  bits  which  are  readily 
available. 

The  important  conclusion  from  this  investigation  is  that 
the  optimum  regime  of  operation  is  around  1  to  1.2  mm 
diameter.   Smaller  holes  will  yield  little  increase  in  gen- 
erated pressure, 
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Tumps  used  in  static  tests. 
Figure  TI-s 


The  explanation  of  this  break  is  not  certain.   It  is 
probably  the  result  of  an  internal  circulation  of  the  fluid 
within  the  larger  diameter  pump  which  results  in  lower 
oressure.   Sea  figure  TI-5. 


Small  Diameter  Pump 


Large  Diameter  Pump 


Figure  TI-5 
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Conclusions 

The  follov.'ing  conclusions  are  drawn  about  the  static 

i 

pressure  produced    by  pumps  with   the  proposed   geometry   of 
chapter   T. 

(1)  The   static  pressure    is   related    to   the   applied 
voltage   by  fi> ~  C(    (~x  J    wiiere   c  depends   upon 
pump  dimensions   and   freon  purity 8      It    is   usually 
in  the    range    0.8   to   1.1    . 

(2)  The    static  pressure   at    constant   voltage    takes   a 
sharp  droo   for  pump   diameters   greater   than  1.6  mm 
if  pump  length   is   held   constant    at    1/16   inch 
(1.59   mm).      It    is   approximately   constant    for 
diameters    less    than  1.6   mm  and   decreases    slowly 
after   the   large   drop    for   larger  diameters.. 
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CHAPTER  ITT 
DYNAMIC  PUMPING 
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DYNAMIC  PUMPING 

If  the  ion  drag  device  Is  to  be  used  as  a  pump,  the 
fluid  must  move,  and  the  ma theme t leal  model  must  take  this 
velocity  into  account.   In  this  chapter  the  theory  of  the 
last  chapter  is  altered  to  account  for  fluid  motion.   The 
generated  pressure  Is  determined  as  a  function  of  fluid, 
velocity.   Then  the  hydrodynamics  of  the  pump  are  investi- 
gated and  the  models  combined  to  give  a  prediction  of  total 
pressure  vs.  velocity.   Measured  pressure-velocity  char- 
acteristics are  presented  for  comparison  and  it  Is  seen 
that  the  theoretical  model  tends  to  overestimate  the  un- 
stable nature  of  the  characteristic. 
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One  Dimensional  Dynamic  Model 


The  assumptions  of  the  model  for  dynamic  conditions! 
are  similar  to  those  for  the  static  model:   charge  relaxa- 
tion time  of  fluid  is  greater  than  pump  transient  time; 
liquid  appears  as  a  conduction;  one  dimensional  dependence 
of  field  and  charge  density;  no  influence  of  lateral  bound - 
ries. 

The  model  is  shown  in  figure  II 1-2  and  the  equations 
to  be  solved  are  given  in  (1)  through  (*0  below.   Equation 
(^ )  has  been  altered  to  account  for  eonvectlve  transport  of 


charge.   Charge  diffusion  is  again  neglected 


E  *  -  V  V 


x  dependence 


x-o 


f 


X, 


E 


dx   =■     £ 


(1) 
(2) 
(3) 


v/ 


Solving   eouatlons    (1)    through    {k )    for   £~<fxj    and    P  Ck)   yield 
(5)    and    (6)    below. 


-  6 


t 


integrate    X-  o    to  X ,     E=  t0    to    E~ 

solve   for    t   Ctc) 

1T£ 


-h 


Ut-t\ 


E  Lx)  = 


**±-i-(F„+y* 


f-  *& 


b 


(**A) 


(5) 
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Poo- 


9A 


[&<<.+  $]* 


6) 


To  connect  charge  density,  electric  field,  and  fluid  pres- 
sure 

d, 


T* 


d%   ~    6 


I  'h 


[& +(<?.+*)*] 


integrate  **  o        ,      *   and     r>  =  0     ,  ^, 


+/f^7*-f.*-£ 


I 


J 


(7) 


Tc  introduce  voltage,  integrate  equation  (5)  using  equation  (2) 


d*     ~    L    ~ 


V 
I 


*♦*.#}  ¥r  +  fr-f)J    -(^f)3-7*, 


(The  ( Fo  +  -j-)  term  was  introduced  as  part  of  an 
arbitrary  constant. ) 

The  voltage  difference  -  U  +  Uo   will  be  written  as 
U-U* v.here     U*~   represents  the  portion  of  the  driving  volt- 
age necessary  to  generate  the  ions  and  \J     is  the  total  volt- 
age across  the  pump. 


(8) 


To  solve  for  pressure  explicitly  In  terms  of  voltege  it 

I 
is  again  necessary  to  make  assumptions  as  to  the  value  of  the 

electric  field  (  or  charge  density)  at  some  point  in  the 
channel.   In  chapter  IT  experimental  results  tended  to  sup- 
port the  space  charge  limited  model,  i.e.  the  assumption 
that  the  electric  field  is  zero  at  the  emitter  in  the  static 
condition.   This  assumption  again  seems  reasonable  if  the 
fluid  is  flowing  in  the  positive  direction,  i.e.  in  the 
direction  in  which  the  ions  tend  to  drag  it.   However  if  the 
fluid  motion  opposes  the  ion  direction  this  assumption  is  not 
reasonable  especially  if  the  fluid  velocity  is  on  the  order 
of  b E  ,  the  ion  velocity  relative  to  the  fluid.   Stuetzer 'sVA  ' 
assumption  was  that  the  current  density  must  be  zero  at  the 
emitter  for  negative  fluid  velocity.   Pickard'^'  introduces 
a  parameter  which  indicates  the  value  of  the  electric  fitld 
at  thr  emitter  and  includes  this  parameter  in  the  pressure 
equation.   However  the  dependence  of  this  parameter  upon 
fluid  velocity  is  not  known,  and  the  assumption  of  Stuetzer 
represents  one  specific  value  of  Plckard's  parameter. 

The  following  development  will  make  use  of  assumptions 
similar  to  those  of  Stuetzer'2'.   The:  results  will  later  be 
compared  with  experimental  results  for  pressure  vs.  velocity 
measurements,  an  investigation  which  was  not  included  in 
Stuetzer' s  reports.   In  light  of  the  experimental  results, 
the  adquaey  of  this  model  for  explaining  pressure  Instabil- 
ity will  be  discussed. 
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Emitter   Field    Assumptions 


£.-  * 


E . 


From  equation   (7)    for       \j  >    0 


f, 


u 


-t- 


■H  \TT  *(tJ   +  *l-hr 


6  6    +  u;  J    -  h 


from  equation  (^A) 


fr-f^f^-  qi^r 


fi^  Y^\ap  va  -*-  Ai 


■w- 


'  ib     Atj        L     a  I  eh 
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-h 


*1  '/ 


1f  i 
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AX     {    6     J 
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(1£l)$  + 


J.   V 


2V 


V  ?    o 


For     v^  0     from  equation   (7) 


A%1,       \/2 


((■ 


»  -  T  (#/  -  ^.    o 


(9A) 


3^ 


i^Y-ii^fT^. 


r 


MiM 
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To  obtain  a  relationship  "between  pressure  and  voltage 
start  with  equation  (8)  written  as  follows: 
for  V  ^  o 


Ea  =    * 


U-  0*  \*9*<    aTfxX'W 


**, 


jt  +iv  r\  (tJ-n 


from    (7) 


>\'A  +  ± 


+ 


it* 


Hf 


a 


5  %/ 


&)*+*]'-&)' 


J 


(i^  [(££-)*,    vj 


3  v 

7T 


-  M 


^  y.    3  / 


-f-    3 


jpEjI  _ 


c-py*/j<if. 
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3  V 

77 


(y-^ 


A£v/ 


^  *  / 


(£&) 


■^  J 


V 


\/^  0 


(10A) 


for     v/^  0» 
eq.    (S) 


E    -   -- 

U-  0*      _ 


**, 


.=2.  J/. 
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from  (?) 


j^f  -  (¥  *  ay  y>  i 


v  ^    o 


A  -- 


^ini 


v  \j  A_y_ 


(lOB) 


V^sD 


e   k*  (      jl 


V 


(11) 

In  order  to  visualize  these  results,  it  is  helpful  to 
normalize  the  variables  as  follows: 

(a)   Pressure  v.  ill  be  normalized  in  terms  of  the  static 
pressures 


-  r> 


~n 


f.-f  (rgj 


? 


(b)  Electric  field  is  normalized  by  defining  -f^  -    ^ 


-^e 


where  £   is  the  average  field  which  would  produce  the  static 
pressure  according  to  the  above  equation 

(c)  Velocity  is  normalized  to  i>  C  ,  the  velocity  oi 
ions  in  static  fluid  under  the;  influence  of  the  normalizing 
field 


V*-- 


x 


(d)   Distance  is  normalized  in  terms  of  the  emitter 
collector  spacing  X  , 
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X 


X 


other  normalizing  factors  follow 

(e)   Current  is  normalized  to  the  static  current 


(f)   Charge  density  Is  normalized  to  the  charge  density 

at  the  collector  under  static  conditions 

& 


R=  n,t 


e  -    x- 


e 


All.  normalized  (dimensionless)  variables  are  vvrltten  with  < 
subscript  ~n   . 

The  equations  for  pressure,  current  density,  electric 
field,  and  charge  density  can  be  normalized  and  are  sum- 
marized as  follov.s: 
Normalized  Variables 


pressure  ^p  7t,  =  ^-—  ,  ^5=  static  pressure 


<  Ci?*) 


distance 


z 


X  n.  -    "x"     »   ^/=  distance    from  em 


ltter   to   collector 


E 


E.    field    t  n  -   "g"       ,       £  =   /fj^ZT 


(15A) 
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Velocity  Vyt=y-,  tt  =  L    £ 


Current  density  a  -ru-    -%=     ,       j-_   i>  ^o^- 
Charge  density   (°t?.  =  "^~  ,   £>  -  _£_   £" 


5  *, 


Normalized   Equations 


o  r»  — 


V    7Z    ^      O 


L(  -f>  7t  +    \/-n)  ^  +     V  tzJ  v  -n  <^    e      (1?) 


E  n    =  7t    Xn    -h      V  tl*    J         -      V-n  V tl  >   0 


(  Qtl    X-n.)        —    V-n  V-n^  o 


(18) 


<° 


72. 


'-¥k)^  Vn  <  a  (15) 


Figures   III-2  through   ITT-5   shov<   graphical   representations 
of  equations    (16)    through    (19). 

Figure    TT1-?   shows    that   the    pressure    increases   to    — — 
of   the    static   pressure    for   large   positive   velocity,    and 
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fells  to  zero  at  v=  -  j,   />  £ .   The  current  density  shown  In 
figure  TII-3  goes  to  zero  at  \s-  -  3  b  £      (which  was  the  ori- 
ginal boundry  condition)  end  becomes  linearly  dependent  upon 
velocity  for  large  V    where  convection  of  charge  is  the 
dominant  transport  mechanism. 

Electric  field  and  charge  density  are  plotted  8s  func- 
tions of  position  along  the  channel  with  the  fluid  velocity 
as  a  parameter.   Tt  is  seen  that  for  V  >  o    the  charge  is 
swept  downstream  toward  the  collector.   This  produces  an 
increased  electric  field  in  the  vicinity  of  the  collector 
since  we  are  operating  at  constant  voltage,  and  the  Integral 
of  the  electric  field  must  remain  constant. 
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Figure   III-2 


Pressure  vs. 
Velocity 
(normalized ) 


v      (normalized   velocity) 
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Figure   III-3 
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Figure  III-4 

Electric  Field  vs.  Position 
(normalized ) 


x    (normalized  do sit  ion) 
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Figure  IIT-5 

Charge  Density  vs.  Position 
(normalized) 


.H  .5     .6,     .7 


x    (normalized  -position) 
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Pressure  difference  observed  across  a  pump  with  fluid 
flowing  is  due  to  the  generated  pressure  used  minus  the 
hy&rodynamlc  pressure  loss  due  to  creation  of  turbulence \ 
beyond  the  oump.   The  pump  can  be  modeled  as  en  orifice  and 
analyzed  using  momentum  conservation. 

To  develop  the  pressure  drop  across  an  orifice  we  assume 

o    1  a. 

that  all  energy  dissipation  occurs 

in  the  region  downstream  of  the 

- 
orjfice,  i.e.  betv:een  sections  1 

and  2.   Applying  the  momentum  Figure  TII-6 

equation  to  the  dashed  control  volume  of  figure  IIJ-6. 
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=    (■?,  -  fp*)   AA"  Cq(vA'V,) 
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From  continiuty 

From  Bernoulli's  equation  the  heaci  loss  is 
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3 


?3 


■3 


Substituting 


H 


L  ,-1 


*  3 


I  - 


A  A 


(20) 


J 


The  flow  contracts  in  flov.ing  through  the  orifice  so 

that  the  effective  area  at  section  c  is  related  to  A,  by 
/,c  =  eA, 

Penlaring  A,  by  cA,  and  using  continiuty 


W 


V,  c  A  =     vo  A 


or   since    Q  -    V&  Aa     (volume   flow    rote) 

a  J-        _£JL     /  Aa 

Hl  -   *3      a/    ITT 


The   vslue   of    c   Is   given  experinentally  as    (8) 

3 


c  =    . «  *  .  3  ip  64^-) 


Tutting  in  dimensionless  form  similar  to  equation  (16) 

A*  J 


*@L    =     T 


A^vt=     C        X         J    )LH where  O  X  (  ''  -£- L ) 

-IPs  A  *  J 


A  tf-n=     b   v^  D=  C        (JJO  (21) 

The  final  equation  for  the  pressure  difference  across 
a  pump  is  the  combination  of  equation  (21)  and  (16).   The 
pressure  difference  of  equation  (21)  will  always  be  a  loss, 

will  always  be  a  pressure  drop  in  the  direction  that  the 
fluid  is  flowing.   Equation  (16)  does  not  give  a  value  for 
the  pressure  when  the  normalized  velocity  is  less  than  -  j  , 
i.e.  when  the  fluid  is  flowing  backwards  faster  than  the 
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ions   are    being  dragged    forward    relative    to   the    fluid.      The 
assumption  will  be   made   that    in  this    case    the   generated 
pressure    is   zero.      The    ions   are   being  washed   back   into   the 
resevoir  behind  the    emitter  and   probably  do   not   act    to   exert 
much  oressure  rise    in  the    forward  direction.      The    resulting 
eauation  is 


-*>■*  = 


'  ^  +  DvS 


/A   "*■ 


Kh.  + 1    -  D 


3     S.    Vvt    *-   O 


V-H    <    -    T 


Vn  >0 


(22) 


A  plot  of  equation  (22)  is  shown  in  figure  JII-7.   The 
slope  of  the  pressure-velocity  curve  is  positive  at  v^-  o . 
This  was  the  basis  for  Stuetzer ' s* 3 )  prediction  of  instability. 
In  the  model  which  he  proposed  the  pressure  was  not  modeled 
for  values  less  than  v^  ^  -— .   If  the  same  equation  as  (22) 

3 

for  x<^^oiE  extended  to  the  region  v-^-i-fL,  the  pressure 
will  decrease  faster  than  the  orifice  pressure  loss  increases, 
and  the  r>umo  would  not  recover  from  the  negative  velocity. 
It  is  necessary  that  the  pressure-velocity  characteristic 
reverse  slope  in  the  negative  velocity  region,  as  well  as 
the  positive  region,  to  have  a  pressure  oscillation  of 
confined  amplitude. 
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Experimentally  jfessured  Pump  Characteristic 

The  experimentally  measured  pressure-velocity  character- 

I 
istic  is  shov.n  in  figure  III-7.   This  is  for  the  purno  whqae 

geometrical  parameters  v;ere  used  in  developing  the  theoreti- 
cal curve  of  figure  III-7. 

It  is  seen  that  the  slope  at  V^  =  o    is  positive  indi- 
cating instability.   The  measurements  were  made  using  an 
experimental  setup  as  shovm  schematically  in  figure  III-8. 
The  flow  rate  was  adjusted  by  changing  the  valve  setting 
and  simultaneous  measurements  of  flow  rate  and  pressure 
difference  across  the  pump  were  made.   Since  the  flow  was 
driven  by  the  height  of  the  resevoir,  the  pump  could  be 
reversed  In  the  enclosure  r-na    similar  measurements  taken 
with  negative  velocity  relative  to  the  normal  pumping 
directions.   A  picture  of  the  experimental  setup  is  shov.n 
in  figure  TTT-9. 

The  measured  characteristic  does  show  the  positive 
slope  at  •=  o  which  is  necessary  for  a  condition  of  unsta- 
ble pressure.   The  pump  for  which  these  measurements  were 
made  oscillated  when  connected  directly  to  8  manometer.  The 
measured  characteristic  rise  for  v >  o      is  less  than  that 
predicted  by  equation  (22)  as  is  the  decrease  in  pressure 
for  v  <  o   . 
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Conclusions 

The  forms  of  the  measured  pump  characteristic  tends  to 
support  the  mathematical  model  used  to  develop  the  theoreti- 
cal characteristic.   The  discrepancy  In  magnitudes  of  the 
humps  just  to  the  left  and  right  of  v  ~  6    can  probably  be 
traced  to  the  inadequacy  of  the  bound ry  condition  applied 
at  the  emitter. 

Although  the  condition  of  zero  electric  field  tends  to 
be  supported  by  the  evidence  presented  in  chapter  II,  this 
assumption  does  not  appear  to  hold  when  the  fluid  has  a 
velocity  greater  than  zero.   If  the  field  at  the  emitter 
became  non  zero,  i.e.  the  charge  density  at  the  emitter  is 
less  than  infinite,  the  generated  pressure  would  fall  as 
shown  by  Pickard(5).   The  generated  pressure  would  not  rise 
to  a  value  ~q~    times  the  generated  pressure  at  zero  velocity, 
and  a  smaller  pressure  rise  than  predicted  by  equation  {??) 
for  small  positive  velocity  would  be  exhibited. 

For  negative  velocity  it  has  been  assumed  that  the 
emitter-  electric  field  asjusts  itself  such  that  no  current 
passes  the  emitter  in  the  negative  direction,  i.e.  that  the 
space  charge  in  the  vicinity  of   the  emitter  continues  to 
build  up  to  produce  the  required  emitter  field.   If  in  fact 
the  space  charge  around  the  emitter  is  not  this  large,  the 
generated  fluid  pressure  would  decrease  to  zero  sooner  than 
predicted  by  the  mathematical  model  and  the  pressure  dip  for 
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velocity  slightly  less  than  zero  would  be  less  than  predicted 
by  equation  [22)  „ 

There  is  no  way  to  predict  how  the  emitter  boundry 
condition  should  change  with  velocity.   A  dependence  of 
emitter  electric  field  could  be  deduced  by  working  backwards 
from  the  observed  characteristic;  however,  this  would  be  a 
circular  argument  and  would  add  no  additional  understanding 
to  the  problem. 

Further  evidence  of  the  validity,  If  not  the  exact 
magnitude  of  the  theoretical  characteristic  is  shown  by  the 
observed  pressure  instability.   The  tneoretical  implications 
and  experimental  observations  relative  to  this  problem  are 
presented  in  chapter  TV. 
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CHAPTER  IV 
PRESSURE  INSTABILITY 


5? 


Pressure  Instability 

The  discrepancies  between  measured  pressure  velocity 
characteristics  and  those  predicted  by  the  theory  of  chapter 
TTI  are  discussed  in  relation  to  the  observed  and  predicted 
motion  of  a  pump- manometer  experiment.   Deficiencies  in  the 
model  are  suggested,  that  may  account  for  these  discrepancies, 
and  areas  of  improvement  are  suggested. 


Pump- Manometer  Instability 

If  a  pump  with  the  pressure-velocity  characteristic  of 
equation  (22)  chapter  III  is  connected  to  a  manometer  a 
pressure  oscillation  should  result.   If  the  manometer  and 
connecting  ducts  are  must  larger  than  the  pump  opening,  it 
can  be  assumed  that  all  hydrodynamlc  losses  ocrur  at  the 
pump  orifice  and  have  already  been  included  in  equation  (22) 
The  difference  of  level  of  the  fluid  in  the  manometers  is  w 
and  will  have  a  static  value  dependent  upon  voltage  and  the 
r>ump  diameter  as  presented  in  chapter  II.   The  differential 
equations  describing  the  motion  of  the  fluid  in  the  manome- 
ters are 


a.  = 


dt 


V 


where  w,  is  the  mass  of  the  fluid  in  the  manometer,  A^  is 
the  cross  sectional  area  of  the  manometer,  and  -p      is  the 
pressure  as  given  by  equation  (22)  in  normalized  form.   The 
value  of  Jh      is  a  function  of  the  velocity  through  the  pump  v^, 
which  is  related  to  the  velocity  of  the  manometer  fluid 
level  V      by 


n  rr} 
P 


r  A 


The  static  pressure  fifap-o)   will  balance  Cj^     static  nnd 
the  equation  can  be  written  as  a  variation  about  u      static 

5^ 


If  equation  (1)  Is  divided  through  by  the  velocity  it 
cpn  be  written  in  the  form 


%  =  4*-  (-  esq  +  fsf~(vf)) 


where  Jb        is  given  by  equation  (22)  of  chapter  Til.   The 
function  ^  {y^\   satisfies  the  requirements  for  the  exist- 
ence of  a  limit  cycle  as  the  solution  to  equation  (2),  and 
this  solution  can  be  determined  by  the  use  of  a  graphical 
t  e chni que  of  Li e na r d  *  ° ' . 

E  xpe  r  i  me  nt 


/;  pump-manometer   experiment    as    shown    )n   figure   IV-1 
with    the    dimensions   giver    in  table    1   was    conducted.      Using 
this    data,    equation    (?)    becomes 

c[v  -  1.  <D6"X  +    53.£  P-h 


J*  V 


(3) 


where  v   is  in  Aec       and  X  in  on.   The  characteristic 
53.5  Dn  is  olotted  in  figure  IV- 2  using  equation  (22)  of 
chanter  TTT.  Viscous  effects  and  other  hydrodynamic  resist- 
ances were  found  to  be  negligible  compared  to  the  pump 
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orifice  loss.   Superimposed  upon  the  characteristic  13  s  plot 
of  the  graphically  determined  limit  cycle. 

The  amplitude  of  the  oscillation  predicted  by  this  cycle 
is  approximately  Hi  cm  with  a  period  of  1?.5  seconds. 

A  plot  of  the  one  period  of  the  motion  actually  ob- 
served is  given  in  figure  TV-3.   The  oscillation  amplitude 
is  approximately  3  cm  with  a  period  of  ^.6  seconds. 

The  experimental  pump  characteristic  given  in  chapter 
TIT  figure  IT  1-8  was  measured  for  the  same  pump  sometime 
after  the  motion  of  figure  IV- 3  had  been  recorded.   Elect- 
rical breakdown  in  the  pump  had  by  this  time  degraded  Its 
performance  somewhat.   In  figure  IV-^  the  measured  charact- 
eristic is  plotted  and  the  same  limit  cycle  techniques 
applied  to  determine  the  motion.   The  predicted  amplitude 
by  this  graphical  construction  is  1.5  cm.   The  observed 
oscillations  at  the  time  the  characteristic  was  measured  had 
an  amplitude  of  approximately  1.2  cm. 
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Pump-manometer   oscillation  experiment. 
Figure    IV-1 


Pump  Diameter 

Ms no  me  t  e  r   A rep 

i 

Fluid  Mass 
Static  Head 
Electrode  Spacing 

Mobility 

Pump- ma no meter  dimensions 
Table   1 


.79    mm 

.18    cm2 

.04    kg 

7.60   cm 

1.59    -Tim 


-7 


57 


X   P, 


Pi 

o 

•H 

*—» 

-P 

-P 

•H 

o 

W 

rH 

o 

ft 

p^ 

CD 

• 

rH 

C\i 

(0 

O 

1 

> 

>> 

> 

C 

H 

>» 

-p 

43 

0) 

<H 

«H 

i-i 

O 

R 

2 

O 

°n 

00 

rH 

r-l 

•H 

0 

— ' 

fo 

> 

58 


Figure  TV- 3 
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Time    (experimental) 
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Conclusions 

The  pressure  oscillation  amplitudes  calculated  on  the 
basis  of  the  measured  pump  characteristic  are  very  nearly 
equal  to  the  observed  amplitudes.   This  indicates  that  the 
mechanism  causing  the  instability  is  increase  in  generated 
pressure  for  positive  velocity  as  predicted  by  the  theory 
of  chanter  TIT  and  given  in  equation  (22)  of  chapter  III. 
The  introduction  of  an  additional  load  in  the  manometer 
system  by  clamping  one  o^  the  tubes  to  a  very  small  cross 
section  caused  the  oscillation  to  stop  and  the  pressure  to 
assume  a  steady  value.   The  effect  is  to  introduce  viscous 
loading  which  eventually  causes  the  characteristic  to  assume 
a  negative  slope  at  velocity  equals  zero,  i.e.  a  stable, 
condition. 

It  had  been  suspected  earlier  in  this  investigation  that 
the  mechanism  of  the  instability  was  a  bulk  instability  in 
the  fluid.   Evidence  of  such  a  bulk  instability  in  the  form 
of  changes  in  current  at  critical  voltages  is  given  in 
chapter  V.   The  onset  of  this  apparent  bulk  instability  is 
at  a  voltage  belov;  the'  operat  lng  point  at  which  the  pressure- 
oscillations  re  re  observed. 

The  discrepancies  in  amnlituc'e  of  the  oscillations  from 
those  predicted  by  theory  are  the  result  of  the  smaller  humps 
in  the  measured  pump  characteristic.   The  comments  in  the 
conclusions  section  of  chanter  IV  apply  here  again. 
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CHAPTER  V 
CHARGE  TRANSPORT 
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Charged  Particle  Motion  Within  The  Fluid 

The  experimentally  determined  curves  of  current  vs. 
voltage  squared  shows  a  very  definite  break  in  slope.   This 
suggests  that  at  this  point  there  is  a  change  in  the  me- 
chanism of  charge  transport  or  a  change  in  the  number  of 
charged  particles  being  produced.   The  possibility  of  the 
onset  of  an  instability  in  the  bulk  is  investigated,  by  look- 
ing at  the  ratio  of  electrical  forces  induced  by  the  motion 
to  viscous  dissipating  forces  (electric  Hartmann  number) 
where  the  break  is  observed. 

P'igure  V-l  shows  a  plot  of  current  ( microampers )  and 
static  pressure  head  vs.  voltage  squared.   The  current 
appears  to  be  proportional  to  voltage  squared  in  two  regions 
with  a  much  larger  slope  in  the  region  U    >  120   (  KVj  o   rp^g 
pressure  also  has  a  slight  break  in  slope  at  approximately 
the  same  -point. 

A  possible  explanation  for  tr.ls  break  Is  the  onset  of  a 
bulk  instability  within  the  fluid  occuring  at  the  voltage 
where  the  break  occurs.   The  criterion  for  the  onset  of  such 
a  bulk  instability  is  the  relative  magnitudes  of  viscous 
dissipating  forces  in  the  fluid  and  electrical  forces  induced 
by  virturc  of  the  fluid  motions. 
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From  equation   {22)    chapter   III    (with  normalizing  vari- 
bles   reinserted) 

/  •=  Z> 
v;here  fij>s     is   the    characteristic   electric    field. 

If  we   assume   a   Poiseuille   distribution  of  velocity  with- 
in  the   cycllndrical   channel   the   pressure   drop  across   the 
channel    is 

?  Ms 


f      = 


R 


f 


d_A£_    =    _£_ 


f 


^     " 


L  ^o 


For  instability  the  pressure  rise  due  to  elect,  forces  must 
increase  faster  than  the  viscous  dissipation  rate. 


.±±s 


> 


l(ok> 


L 


>      / 


r—  (1) 

From  equation   (9)    chapter   II   the   value    of     $fs     is   derived 

by   measuring  the    slope   of   the   n     vs.   (J*  curve. 


^C 


For  the  experiment  shown  in  figure  V-l  the  value  of  C  is  0.73 
and  /  ?l£s  =  7.  £5"x  to    •      For  this  experiment  ^  £ 43=  Lp  .   The  value 
of  b  la  determined  by  simultaneously  measuring  current  and 
pressure  vs.  voltage  squared. 


6*» 


is  ,    Q  J_  (  u  )* 

5     L5fUy 


x  = 


ca  *(£)* 


By   solving   for   C   In  the    first   equation,    using  this   value   in 
the   second   equation  and    solving   for  b  a   value    of  l  =  .3of Klt>'7 £^~^if- 
results.       (The   value   given  by    Stuetzer   for   freon  113  was 
7X/o~7  .)      Using  the   value    t  =,  3c/X  /D'1  and 


-   .  1    K    /o 


t  J*** 


/(o  b 


A±S 


L 


=  ;9I 


Since  f^_±s__    represents  an  electric  field  the  express- 

^     €  e  U 

1  on  en  the   left   of    Inequality    (1)    has    the    form    ~ZuH>  or  £n 

electric  Hartraann  number.      Schneider   and  fcfatson* 9   ;   have 
shown  that   the    critical   value    of   this   parameter    for  a   planar 
geometry  v.as   99  •  5   for  the   onset   of   the   bulk   instability.      The 
value    of   191   for   this   experiment    is   the    same   order  of   magni- 
tude  and    subsequent   experiments   have    shown  that    in  fact   the 
critical   voltage    for   the   break    in  the    current-voltage    sc'uared 
curve    is   practically   independent   of  Dump   length  and  transverse 
dimensions   as  would  be    implied  by   the    form    — -jr        . 

For   the    oscillating   manometer  experiment   described   in 
chapter    TV   the    voltage   was   above    the   value   at    i hich   this 
current   break    evidence    of    bulk    instability    occurs.      It    is 
probable   that   the    numn   fluid    is   nearly   always    in  the    turbu- 
lent  condition  produced   by-  a   bulk   instability,    and    that    this 
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type  of  Instability  is  not  the  mechanism  causing  the  pres- 
sure oscillations. 
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CHAPTER  VI 
WICK   DESIGN 
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Pes 3  gn  Parameters 

The  objective  of  investigation  of  the  proposed  ion  drag 

i 
pump  geometry  was  to  develop  enough  information  about  their 

performance  that  an  electrohydrodynamlc  wick  could  be  design- 
ed using  arrays  of  these  pumps.  In  order  to  produce  an  oper- 
able wick  the  design  parameter  which  must  be  supplied  are 

(1)  Electrode  thickness 

(2)  Electrode  spacer  thickness 

(3)  Pump  hole  diameter 

(b)      Pump  spacing  within  an  array 

(5)   Spacing  of  the  arrays  relative  to  each  other 

As   stated  previously,  the  thickness  of  electrodes  has 
little  effect  upon  pump  performance.   For  most  experimental 
studies  brass  with  thickness  .021  inches  was  used.   The  more 
important  influence  uoon  performance  is  the  size  of  burr  pro- 
duced in  drilling  which  depends  upon  electrode  hardness,  bit 
sharpness,  and  rotational  speed  of  the  drill.   A  burr  which 
protrudes  approximately  „ ?mm  from  the  electrode  surface  seems 
to  work  well. 

Electrode  spacing  for  most  experiments  was  achieved  with 
1/16  inch  plexiglass  betv.een  the  electrodes.   This  spacing 
was  found  to  be;  best  with  readily  available  stock  plexiglass 
and  war,  the  result  of  two  considerations: 

(a)   Smaller  spaces  such  as  1/32  inch  plexiglass  often 
resulted  in  electrical  breakdown  at  very  low  vol- 
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tages.   This  probably  because  the  spacing  is  becom- 
ing almost  as  small  as  the  burr  projection.   Since 
the  burr  varies  from  one  pump  to  another  the  pro- 
•  bability  that  a  portion  of  the  burr  may  extend 
across  the  spacer  at  some  point  is  high  and  break- 
down is  quite  likely; 
(b)   Larger  spaces  such  as  1/8  inch  require  that  the 
voltage  be  doubled  to  produce  the  same  generated 
pressure  in  the  pump.   Voltages  result  (25-30KV) 
which  may  cause  breakdown  outside  the  pump,  and 
were  not  as  easily  produced  because  of  available 
power  supplies. 
Pump  hole  diameter  selection  should  be  the  result  of  a 
tradeoff . between  higher  generated  pressures  for  smaller  dia- 
meter rjunros,  and  lower  hydrodyramic  losses  in  the  pump  it- 
self for  larger  diameter  holes.   However  the  experimental 
Investigation  of  generated  static  pressure  vs.  hole  diameter 
presented  in  chapter  ITT  clearly  indicates  thtt  the  optimum 
region  of  operation  is  with  diameter  just  below  1.6mm.   This 
insures  that  the  pump  will  operate  with  the  higher  generated 
pressures  just  above  the  break  in  the  curve. 

Spacing  of  holes  within  an  array  snould  be  as  close  as 
possible  to  deliver  the  maximum  flow  rate  with  the  smallest 
wick,  and  to  reduce  hydrodynamic  losses  as  much  as  possible. 
The  effect  area  seen  outsice  the  orifice  (area  A,  in  equation 
20  chanter  TIT)  becomes  less  if  the  pumps  are  placed  closer 
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together  and  the  orifice  losses  are  consequently  less.   The 
factor  which  determined  how  close  together  the  pumps  could 
be  placed  was  drill  accuracy  w hi ch  required  that  the  center 
to  center  spacing  be  about  two  diameters.   Closer-  spacing 
might  result  in  cracks  in  the  material  between  the  holes 
even  if  more  accurate  positioning  of  the  drills  were  employed. 
No  larger  effect  upon  generated  pressure  was  observed  with 
spacings  which  varied  from  2  diameters  to  ?0   diameters. 

Spacing  of  pump  arrays  relative  to  one  another  is  de- 
pendent urion  the  pressure  head  which  is  maintained  at  the 
flow  rate  at  which  thf  array  is  operated.   For  wick  operation 
it  is  only  necessary  that  one  stage  lift  the  fluid  vertically 
through  the  next  stage  above  it.   The  usual  objective  of 
experiments  was  to  operate  at  a  flow  rate  such  that  approx- 
imately 1.5  centimeters  head  was  developed.   Spacing  of  arrays 
approximately  one  centimeter  apart  insured  that  each  array 
would  be  able  to  deliver  the  fluid  to  the  next  stage.   Very 
close  spacing  will  result  in  some  backwards  drag  of  ions 
and  a  pressure  loss  between  stages.   This  effect  is  minimum 
because  of  the  open  reservoirs  between  arrays.   The  ions 
ere    not  directed  to  give  maximum  drag  effect  as  they  are  in 
the  Dumps.   Tf  a  spacing  of  at  least  %   inch  is  maintained 
between  arrays,  the  pressure  loss  is  negligible. 
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An  Experimental  wick 

A  vilck  of  3  arrays  was  constructed  as  shown  schemati- 
cally In  figure  VT-1.   Pictures  of  the  arrays  and  the 
assembled  wick  are  in  figure  VI-2. 

The  wick  dimensions  were  as  follows: 

Number  arrays  3 

Number  of  holes  per  array  25 

Hole  diameter  1.40  ram. 

Center  to  center  hole  spacing      J.O   ram. 
Array  electrode  spacing  l/lo   inch 

Electrode  thickness  .021   inch 

Spacing  between  arrays  3«5  cm. 

The  spacing  between  arrays  3»5  centimeters  to  per- 
mit observation  of  the  flow  and  possible  electrical  break- 
down if  it  took  place  in  any  of  the  pumps.   A  spacing  of  1 
centimeter  would  have  been  more  than  adequate  to  isolate  the 
electrodes  of  adjacent  arrays. 

Performance 


The  array  was  tested  using  freon  doped  with  copper 
nanhenate  at  two  voltages.   At  10  kilovolts  the  static 
oresnure  developed  was  5.9  centimeters  head.   With  the  load 
valve  opened  such  that  the  flow  rate  was  2.1  cra-^/sec,  the 
pressure  was  ^4.5  em  head.  TYie    flow  rate  of  2.1  criK/sec  was 
the  maximum  reading  on  the  flov.  meter  used;  however  by  de- 
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creasing  the  load  and  pressure  across  the  load  larger  flow 
rates  were  possible.   At  12.5  kilovolts  the  pressure  at  zero 
flow  rate  was  2k . 6  centimeters  and  the  pressure  at  2.1 
cm^/sec  flow  rate  was  22. 8cm. 


Evaluation 

Based  upon  the  above  wick  experiment  along  with  several 
other  similar  experiments  the  following  conclusions  are 
drawn: 

(1)  Large  flow  rates  can  be  achieved  by  building  wick 
arrays  with  many  holes  in  parallel. 

(2)  The  pressure  generated  across  the  arrays  is  more 
than  sufficient  to  raise  the  fluid  to  the  next 
array  above  it.   Enough  arrays  placed  in  series 
hydrodynaraically  (and  in  parallel  electrically) 
can  raise  the  fluid  to  any  arbitrary  height. 

(3)  Pumps  have  been  operated  for  several  hours  with- 
out noticeable  deterioration  if  the  voltage  is  kept 
low  enough  to  prevent  electrical  breakdown. 

(J|)   The  pump  configuration  proposed  in  chapter  I  is 

extremely  simple  to  fabricate.   Satisfactory  pump 
arrays  can  be  produced  without  special  care  and 
with  a  variety  of  materials. 
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Three   array  v;ick   experiment 


Figure   VI-1 
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'hree   array  wick    experimental    setup. 
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.Single  array  fr on  wick  experiment. 
Figure  VI-2 
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RECOMMENDATIONS  FOR 
FURTHER  INVESTIGATION 
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Pe commendations  For  Further  Invest igat  1  en 

Improvement  of  the  static  ma the ma tics 1  model  might  be 
investigated  by  trying  to  determine  the  field  pattern  or  a 
better  picture  of  the  charge  particle  trajectories  in  the 
Dump.   This  might  be  done  using  an  optical  technique  al- 
though the  size  of  the  pumps  makes  this  difficult.   A  set  of 
current  measurements  at  different  pump  diameters  at  constant 
■oump  length  and  vnltage  might  reveal  more  about  the  oossibil- 
ity  of  an  internal  fluid  circulation  in  large  diameter  pumps 
as  suggested  in  chapter  II. 

More  accurate-velocity  measurements  at  very  small 
velocity  (both  positive  and  negative)  should  be  made  and  a 
comparison  of  motion  predicted  by  this  characteristic  to 
recorded  motion  should  be  made.   Further  study  of  corona 
emitters  in  moving  liquids  might  indicate  a  possible  refine- 
ment of  the  emitter  field  assumption  end  improve  the  pre- 
diction of  the  observed  experimental  characteristic. 

Other  fluids  such  as  transformer  oil  should  be  tested 
further  in  the  proposed  pumps.   Some  experiments  by  the 
author  shoved  that  the  pumps  worked  using  transformer  oil, 
but  extensive  expc  riments  vers  not  carried  out  because  of 
dlff icu"1  tity  in  cleaning  the  apparatus  for  further  freon 
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experiments  after  using  oil.   The  Possibility  of  using  these 
pump  arrays  to  circulate  oil  in  transformers  provides 
motivation  for  such  investigators. 
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APPENDIX  I 
FABRICATION  TECHNIQUE 
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Fabrication  Technique 

The  objective  of  developing  a  simple  geometry  which) is 

i 
easily  fabricated  is  fulfilled  by  the  geometry  proposed  in 

chapter  I.   To  Insure  that  the  technique  used  in  fabricating 
the  pumps  used  in  experiments  reported  in  this  thesis  may  be 
reproduced  by  other  Investigators  or  used  in  constructing 
wicks  for  application,  the  procedure  used  is  described  here. 

The  factor  which  affects  performance  of  the  pump  most 
is  the  size  of  the  burr  which  is  produced  in  the  drilling 
process.   This  depends  upon  the  hardness  of  the  metal  plate 
electrodes  used,  the  sharpness  of  the  drill  bits,  and  the 
rotational  speed  of  the  drill.   Thickness  of  the  metal  plates 
seems  to  have  only  a  very  slight  effect.   Satisfactory  pumps 
were  produced  with  several  different  brass  electrodes  and 
also  with  aluminum.   The  thickness  of  the  insulating  elect- 
rode determines  the  magnitude  of  the  electric  field  at  a 
given  voltage  and  therefore  the  pressure  generated.   Exper- 
ience with  many  experiments  showed  that  a  spacer  less  than 
1/16  Inch  thick  tended  to  arc  at  lev:  voltages  in  many  pumps 
because  of  the  irregularities  of  burrs  produced.   Hence  1/16 
inch  insulating  spacers  were  adopted  for  almost  all  subquent 
experiments.   Use  of  this  thickness  produced  satisfactory 
pumps.   For  a  spacer  thicKness  of  1/1 b  inch  a  burr  height  of 
approximately  0.?  mm.  was  found  to  work  well.   This  may  be 
estimated  by  inspecting  the  size  of  the  burr  on  the  electrode 
from  which  the  bit  emerges.   The  burr  produced  betv.een  the 
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plates  is  obviously  the  one  which  is  desired  but  inspection 
of  the  other  one  will  give  indication  of  size.   A  very  large 

burr  produces  by  a  dull  bit  may  cause  the  metal  electrodes 

I 
to  be  too  close  together  inside  the  hole  and  produce  elect- 
rical breakdown  at  low  voltage.   A  very  small  burr  will  not 
produce  charge  as  well  and  pumping  is  limited.   However  it 
should  again  be  stated  that  satisfactory  performance  can  be 
achieved  without  an  optimum  size  burr. 

Electrode  olates  were  secured  to  the  insulating  spacer- 
be  fore  holes  were  drilled.   Attempts  to  drill  the  holes 
separately  in  the  electrodes  and  spacer  with  allignment 
afterward  were  seldom  successful  and  required  considerably 
more  effort  than  securing  the  electrodes  first.  The   plates 
and  plexiglass  spacers  were  roughened  with  sandpaper  and 
washed,  and  the  electrodes  were  glued  using  epoxy.   Typical 
pumps  with  single  holes  for  static  pressure  experiments  are 
shown  in  figure  II-*J  along  with  the  enclosure  in  which  they 
were  placed  to  contain  the  fluid.   Electric;- 1  contact  was 
made  with  the  plates  by  spring  electrodes  in  the  enclosure 
vhich  v ere  connected  to  an  external  high  voltage  power  supply. 
The  insulating  spacer  in  the  pump  sandwiches  was  extended 
approximately  h   inch  outside  the  metal  electrodes  on  all 
sides  to  prevent  arcing  outside  the  pump.   Wick  arrays  with 
many  pumps  are  shown  in  figure  VT-2. along  with  the  similar 
enclosure  in  vhich  they  were  tested. 
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